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Nuclear quantum transport for barrier problems 

Christian RuMMEL and Helmut HOFMANN 
Physik- Department der TU Miinchen, D-85747 Garching, Germany 

A method is presented which allows one to introduce collective coordinates self- 
I consistently, in distinction to the Caldeira-Leggett model. It is demonstrated how the par- 

, tition function Z for the total nuclear system can be calculated to deduce information both 

■ on its level density as well as on the decay rate of unstable modes. For the evaluation of Z 
' different approximations are discussed. A recently developed variational approach turns out 

superior to the conventional methods that include quantum effects on the level of local RPA. 
Dissipation is taken into account by applying energy smearing, simulating in this way the 
coupling to more complicated states. In principle, such a coupling must depend on temper- 
' ature. Previous calculations along another microscopic approach show this fact to imply an 

intriguing variation of the transport coefficients of collective motion with T. The relevance 
of this feature is demonstrated for the thermal fission rate and for the formation probability 
of super-heavy elements. 

T— I ■ 

> : 

O . §1. Problems of the Caldeira-Leggett model 

I from nuclear physics point of view 

. The two most prominent examples of nuclear physics where processes are gov- 

I erned by (iso-scalar) motion across potential barriers are fusion and fission. In both 

■ cases dissipation may play an important role. For large excitation energies the pro- 
. cess is dominated by thermal activation. With decreasing energies quantum effects 

become more and more important until at zero temperature mere quantum tunnelling 
is the only contribution at sub-barrier energies. 
ip3 . Dissipation in the collective degrees of freedom describes flow of energy to a large 

^ I set of other, fastly relaxing degrees of freedom, commonly referred to as the "heat 

• bath" although in nuclear physics this notion requires greatest care. Sometime ago 

. Caldeira-Leggett have suggested a modeP which allows for an exact treatment of 

a I the bath degrees of freedom. This becomes possible at the expense of an oversimplifi- 

cation of both the bath itself as well as of its coupling to collective motion. Whereas 
this assumption does not seem to be very restrictive for applications in condensed 
matter physics, it violates essential requirements for nuclear collective motion.'^ Let 
us just mention the most serious deficiencies. 

• In the nuclear case the collective degrees of freedom (CDOF) are not indepen- 
dent of the nucleonic ones, viz of the heat bath. Already the simplest condition 
on self-consistency implies that the shape of the mean field in which the nucleons 
move must vary in non-linear way with the CDOF. 

• Even in the simple case of small amplitude collective motion, the frequencies are 
determined by RPA-like secular equations which differ essentially from those of 
the Caldeira-Leggett model.!^ 

• Related to this issue is the fact that the Caldeira-Leggett model is unable to 
make predictions for the evaluation of the transport coefficients for collective 
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motion like potential energy, inertia and even for the strength of friction. All 
of them depend sensitively on both the collective coordinate as well as on tem- 
perature. 

Without invoking the Caldeira-Leggett model so far quantum effects in nuclear 
dissipative dynamics were mainly treated applying real time propagation,^''^' with 
the exception Whereas the former approach is based on the deformed shell 

model in the novel approach one starts with a typical nuclear Hamiltonian of two 
body nature before CDOF are introduced, as shall be explained below. 

§2. Self-consistent dynamics for separable two body Hamiltonians 

Suppose we are given the Hamiltonian with separable two body interaction of 

k 

n = H + -FF with A;<0. (2-1) 

It may be understood as one term in an expansion of the two body interaction into 
separable terms. For transport problems this ansatz has been used before apply- 
ing the Bohm-Pines method of introducing collective coordinates, se^ with further 
references given therein. These applications were based on a real time approach, im- 
plying that this method is limited to temperatures above a certain Tc. Here, we want 
to follow an imaginary time approach .'^'^^'EJ Different to the Bohm-Pines method 
the collective variable is introduced through a Hubbard-Stratonovich transformation 
in the functional integral for the partition function of H2-l() we arrive at the one body 
Hamiltonian 

Hus = H + q{T)F + ^^q\T) (2-2) 

where the collective coordinate q depends on the imaginary time r = . . .h/T and 
has been introduced self-consistently. After expanding the periodic collective "path" 
into a Fourier series g(r) = qo + X^^^o 9^ exp(ifj.'7") where Ur = (27r//i)Tr are the 
Matsubara frequencies the final form of the functional integral for the partition 
function reads 

Z{T) = I dqo exp(-.FSPA(T, q,)/T) C{T, qo) . (2-3) 

The go-integral represents the thermal fluctuations, while the quantum fluctuations 
with amplitudes qr that become important at low temperatures enter the remaining 
path integral for the "dynamical" corrections 

C{T,qo)= y"p', exp(-^i|^) (2-4) 

where the "action" can be expanded like 

A{T, qo) = Xrs{T, qo)qrqs + PrstiT, qo)qrqsqt + cfrstuiT, qo)qrqsqtqu + ^(^r) • (2-5) 

The exponential in H2-3() represents the partition function of the static part of the 
Hamiltonian H2-2() with J-^^^iqo) being the corresponding free energy. The simplest 



Quantum transport over barriers 



3 



approximation to the dynamical corrections consists in neglecting them entirely: 
C^^^{T,qo) = 1 represents the classical limit and has been given the name Static 
Path Approximation (SPA) in the pastP 

The coefficient A^.^ of the second order part of (|2-5() can be expressed in terms 
of the response function x{T,qo,uj) that describes the changes of the average (F)^ 
for small variations of the path 

Ks{T, qo) = (1 + kx{T, qo, uj)) dr-s = Yi T^" — 2 i \\ ^^-s (2-6) 

Here, the o;^ are the local RPA frequencies and hjJu = €k — ei are the energies of 
the intrinsic excitations associated to the static part of the Hamiltonian H2-2() . In 
cases where Xrs is positive and sufficiently large a truncation of (|2-5|) after the second 
order will do and the remaining gj,^o"™tegrals in D'q can be evaluated exactly. One 
ends up with the result of the Perturbed Static Path Approximation (PSPA) 

C''^^HT,qo) = l{X;HT,qo) (2-7) 

r>0 

which takes into account quantum fluctuations on the level of local RPAI^'^^'EJ 

As Xrs depends on temperature the condition Ai > for the breakdown of this 
approximation due to diverging integrals can be transformed into the requirement 
T > Tq, where for the "crossover temperature" (note the similarity to the case of 
dissipative tunneling'^) Tq < Tc holds true."^ 

The coefficients p^st and Urstu of H2-5() have been calculated irP where the break- 
down of the approximation has been shifted to Tq/2 by treating the qi- and q2- 
integrals in Vq to fourth order analytically. Meanwhile an approximation has been 
developed^ that takes over basic ideas of a variational principle suggested for one- 
dimensional systems and.^ Choosing a suitable Af^ the correction factor H2-4|) 
can be rewritten as the average defined with respect to the weight exp(—Af^/\k\T) 

C{T,qo) = (^exp — 1 ) . (2-8) 

In order to determine the best approximation to the equilibrium quantity 2{T) 
the convexity of the exponential (e~^) > e~^^'^ can be used such that the following 
inequality can be exploited: 

C{T, qo) > C^-{T, qo) = exp ("^^^^^f^^) (2-9) 

For Af^ the PSPA form is chosen but with the RPA frequencies uj^{T, qo) replaced by 
a number of adjustable parameters f^^{T, qo). Maximizing the right hand side of (|2-9j) 
with respect to i7^(r, go) gives an optimal approximation to the partition function 
(|2-3|) . In the left panel of fig. Q we demonstrate at the example of the free energy of the 
exactly solvable Lipkin-Meshkov-Glick model (LMG) that the variational approach 
gives excellent results and is applicable even at very low temperatures. The classical 
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Fig. 1. Left: Free energy of the LMG in various approximations. The variational principle gives 
very good results up to /3 = 1/T = 5MeV~^ whereas the PSPA brakes down at Po = l/To = 
1.13 MeV~^. Right: /3-dependence of the quantities + <-l>rpa{<1o ~ 0) and + f2^^^{qo — 0). 



SPA is good only at high temperatures. Inclusion of quantum effects on the RPA 
level via PSPA improves the approximation at high temperatures considerably but 
breaks down at Tq. The approximation ofj^ called ePSPA, behaves well in the 
crossover region T « Tq. In order to illustrate the reason for the applicability of 
the variational principle at low temperatures we plot in the right panel of fig. ^ the 
/3-dependence of the quantities I'l + a;^p^(T, go = 0) and + f2^^^{T,qQ = 0) that 
determine the sign of Ai and its variational analog via H2-6|) . The former quantity 
becomes negative at /3o = l/?o implying a breakdown of the PSPA, whereas the 
second one stays positive. 

§3. Introduction of dissipation and microscopic transport coefficients 

So far dissipation has not been taken into account. The coefficient Ar(T, go) 
of (|2-6() contains the response function = x'('^) + ^x"i^)- Iii the independent 
particle model the dissipative part typically has the form (for tj > 0) 

x" {^) = fki - ^ki) (3-1) 

k>l 

with strengths fki- We know that (at least at not too small excitations) the simple 
states of the independent particle model are coupled to more complicated ones via 
the residual interactions. The spectrum becomes more dense and the self-energy of 
the individual states acquires a finite width such that one ends up with a continuous 
function x"(^)- We do not go into these details^' here but for simplicity mimic such 
effects by smearing out the spectrum by hand.^' To such an averaged spectrum we 
fit the Lorentzian 

x"{T,qo,u;) = F\T,q,) (- w^,^ " ^ -^)] (3-2) 

\{uj-£{T,qo)Y + [r{T,qQ)/2Y J 
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and extract the three parameters -F^, £ and F. In this way the full response function 
may be interpreted locally as the one of a damped harmonic oscillator: 

X{T, qo,lmz / 0) = Xosc = , ..^ ^ ^ 2 frp ^ "FT^F ^ 2 (^'^^ 

The transport coefficients for the local inertia M, the local nucleonic frequency Wnuci 
and the local damping F can easily be related to the fit parameterJ^ and depend on 
temperature T and the collective coordinate qq simply because the original response 
function has those dependencies. One can easily convince oneself that the solution 
of the secular equation 1 + kxii^) = for the collective frequencies now has a finite 
imaginary part, indicating that collective motion is damped. The big advantage 
of the method suggested here is that all three transport coefficients M, Wnuci and 
r are taken from the same microscopic theory, instead of piecing together various 
macroscopic pictures. 

In order to illustrate the drastic differences between our microscopic approach 
and the frequently used combination of wall friction, liquid drop stiffness and ir- 
rotational flow inertia we plot in flg. |21 the temperature dependence of some ratios 
of the transport coefficients for the barrier of ^^^Th. In the macroscopic picture 




1 2 3 1 2 3 
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Fig. 2. Comparison of macroscopic (dashed) and microscopic (symbols) ratios of transport co- 
efficients for the barrier of ^^''Th. The fully drawn line shows the simple analytical formulas 
proposed in^ Left: F = y/M = ft/rkin, often called (3 in the literature. Right: The effective 
damping 77. The dotted line indicates the change from under-damped (77 < 1) to over-damped 
(77 > 1) motion. 



the ratio F = 7/M does not depend on temperature and is much larger than in 
the microscopic picture where it increases with temperature and vanishes at T = 0. 
Within this macroscopic picture the effective damping strength r] = 7/(2y^|C|M) is 
strongly over-damped implying the applicability of the Smoluchowsky limit. For the 
microscopic transport coefficients, on the other hand, we observe a transition from 
under-damped motion (t/ < 1) to over-damped motion (r/ > 1); note that in the for- 
mer case the inertia plays an important role. Even at T = 3 MeV the Smoluchowsky 
limit is not fully reached. 
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§4. Application I: Thermal fission 

Applying a saddle point approximation to the go-integral in l\2-'6^ we can calculate 
the rate of thermal fission from the imaginary part of the free energy from the 
formulfP*^ 

R{T) = -^^lmT{T) for T > Tq . (4-1) 

In situations where the potential is given by just one minimum located at q = qa 
and just one parabolic barrier at q = qb we obtain 



'^-^iW^l (v^S - e---""- . , ,,2) 

In order to simplify notation we have omitted the T-dependence everywhere. Besides 
the Arrhenius factor there appear non-trivial ones which are worth to be discussed in 
detail. The first two factors can be rewritten in the form zu{qa)/2-K ■ M (qa) / M (qb) 
where the attempt frequency at the minimum and the ratio of the inertias at the 
minimum and the barrier appear. The term in brackets is the famous Kramers 
factoiEHJ 

that takes into account the decrease of the decay rate due to damping. The 
last factor represents the increase of the rate due to quantum corrections. It has been 
calculated for the example of ^^^Th. While^^ use a Caldeira-Leggett type 

model, microscopic temperature- and coordinate dependent transport coefficients are 
used in.lSJ Its contribution turns out to be significant at low temperatures (but with 
T > To). 




'0 0.5 1 1.5 2 2.5 
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Fig. 3. T-dependence of Kramers' factor for macroscopic and microscopic transport coefficients. 

IrP^ the importance of microscopic transport coefficients for the description of 
thermal fission has been examined in more detail. It has been found that the ex- 
perimentally observed phenomenon of "onset of dissipation'^^ cannot be explained 
by the macroscopic set of transport coefficients. In order to illustrate this statement 
we plot in fig. |31the temperature dependence of Kramers' factor. Using macroscopic 
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transport coefficients no "onset of dissipation" can be seen at all, whereas for micro- 
scopic ones it appears around the T ?a 0.5 MeV, exactly where pairing disappears. 

§5. Application II: Formation probabilities in fusion 

Now we like to turn to a more recent application of microscopic T- and go- 
dependent transport coefficients to the probability of compound nucleus formation 
in production of super-heavy elements.^ We do not study the initial phase of the 
reaction explicitly but make the assumption that due to dissipation the sum of the 
incident energy E and the assumed barrier height Eb = 10 MeV is distributed into a 
remaining kinetic energy and a thermal excitation according to 

Ef'' = R{E + Eb) and E^ = {I - R){E + E^) . (5-1) 

For the parameter R we compare the effect of the two assumptions i? = and 
= 0.25 on the final formation probability 

^oo = \ erfc {-^=\ (5-2) 



2 \ y2X;°° 



with 



and 



ioo 



'p: + ./^£l wtae ^_^J_llp± ,5.3) 



2i:,^ = ""''^''l, + (^)' {^M - ^f) . (5-4) 

In fig. 0] we compare the dependence of the formation probability on incident energy 




100 

E in MeV 



200 



Fig. 4. Formation probability as a function of incident energy. The initial conditions Eppito) = 
Epp'^ of thermal equilibrium for the momentum distribution, Eqq{to)Spp{to) = 3ft^/4 and 
Eqpito) = are used. Clearly these fluctuations are compatible with quantum mechanics. 



for microscopic transport coefficients with that of macroscopic ones. In both cases 
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the upper curve corresponds to finite initial kinetic energy {R = 0.25) whereas the 
lower one represents i? = 0. Most striking is the fact that for not too high incident 
energies the weak damping in the case of microscopic transport coefficients results 
in formation probabilities that arc orders of magnitude larger than those of the 
macroscopic picture. In addition to that - as expected for the case of under-damped 
or slightly over-damped motion ry^ ?a 0.5 . . . 2 - the microscopic formation probability 
is extremely sensitive to the initial kinetic energy. Such a dependence is very weak 
in the case of macroscopic transport coefficients. The reason is found in the fact that 
in this case motion is strongly over-damped 77(, ~ 5 . . . 6. 



§6. Summary 



We have developed a method for the description of dissipative collective motion 
in finite nuclei that does not suffer from the deficiencies of the Caldeira-Leggett 
model. This method enables one to calculate microscopically on the same footing 
transport coefficients for inertia, damping and stiffness, including their dependence 
on temperature and collective coordinate. Comparison with the macroscopic picture 
of the liquid drop model and wall friction shows significant differences, resulting in 
much weaker damping at low temperatures in general. These effects are shown to 
be important in thermal fission as well as in the formation of the compound nucleus 
in super-heavy element production. 
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